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Two-photon absorption (TPA) spectra of push-pull chromophores are described in terms of a 
two-state model accounting for electron-vibration coupling. Vibrations have two main effects in 
O^l ' TPA spectra. The most obvious one is the appearance of a vibrational structure in the spectrum; 

in this respect we find large Herzberg- Teller contributions. The second effect was not recognized so 
far: vibrational states contribute a new channel to TPA process, that shows up with a blue-shift 
h-j ' and a distortion of the spectrum. Vibrational-channel contributions to other NLO responses are 

, shortly discussed. 
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g ! I. INTRODUCTION 

The large non-linear optical (NLO) responses that characterize 7r-conjugated materials make them good candidates 
for innovative applications. Among other properties, large two-photon absorption (TPA1 cross-sections are promis- 
ing for different applications, ransdn*; from optical storage to biological optical imagingEI. As a consequence, both 
experimental and theoretical work0~Q have received strong impetus and several interesting strategies have been pro- 
posed to improve NLO responses in generalBtl and TPA cross-section in particulaitluj. A presentlj-iiighly debated 
£>"v topic concerns the importance of vibrational contributions to NLO responses of organic materialsliiTO. The cou- 
| pling of electronic and vibrational degrees of freedom (e-ph coupling) is very important in conjugatedjaolecules and 
i ^~f , polymeraLj and is expected to originate large effects in systems with intrinsically non-linear behaviorEJ3~t3. General 
agreement is presently emerging about the importance of vibronic contributions to static NLO responses in conjugated 
t-H , materialsOOESlia. This can be easily understood: the static applied electric fields affect the electronic properties 
of the chromophore by affecting its 7r-electron distribution. Due to e-ph coupling, the applied fields then also act 
on the vibrational states by displacingjauclei, and this in turn reflects on the electronic distribution, in a feedback 
mechanism that amplifies non-linearity£3. 

The situation is much more complex at optical frequencies. A fairly common statement is that nuclei are too 
massijze to respond at optical frequencies, so that their contribution to NLO responses is irrelevant. In a recent 
paperE3, negligible vibrational contributions to non-resonant NLO responses have been demonstrated. Whereas for 
actual molecules with electronic excitation energies of the order of 1-2 eV and vibrational energies of about 0.1-0.2 eV 
it is difficult to find out a frequency range far enough from both electronic and vibrational resonances as to safely apply 
the approximations in Ref.Ej, one indeed does not expect really large effects of e-ph coupling in truly non-resonant NLO 
responsesli-X As far as resonant processes are concerned, however, the above statement does not apply. Quite obviously 
it leads to erroneous conclusions if applied to incoherent measurements, i.e. to experiments where a 'prepared' state 
■ is allowed to relax for some time before further processing. Just as an example, we have already proved large, and 
non-trivial effects of the coupling .between electronic and slow degrees of freedom in time- resolved fluorescencalj and 
. femtosecond hole-burning spectraE3 of push-pull chromophores. But also in coherent resonant experiments (one- and 
two-photon absorption being simple examples) e-ph coupling is important. Just as in one-photon absorption (OPA) 
spectra, also in non-linear spectra e-ph coupling is responsible for structured, vibronically resolved band-shapes. 
Large effects of e-ph coupling have been recognized within the Condon approximation in TPA spectra of push-pull 
chromophoreaj, and, within the same approximation, TPA, third harmonic generation (THG) and non-degenerate 
four-wave mixing spectra of the conjugated polymer polydiacetylene have been shown to be-non-trivially affected by 
e-ph couplingjij. Recent ab initio calculations of TPA spectra of a few conjugated moleculestrc! showed sizable and, in 
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some cases even large, Herzberg- Teller (HT) corrections, demonstrating once more the non-trivial role of e-ph coupling 
in conjugated systems. 

The effects of e-ph coupling in NLO spectra calculated within the Condon and HT approximations are important 
and have to be included in any reliable spectral simulation. However they are not dissimilar in principle from standard 
effects in linear absorption spectra. In fact, within Condon and HT schemes, e-ph coupling only leads to a partitioning 
of the oscillator strength of a given electronic transition into several vibronic states, or occasionally (within HT) to 
a borrowing of intensity from an allowed transition to a forbidden one. Of course the spectral consequences of 
these well-known phenomena will be more and more complex and important the higher is the order of non-linearity, 
and the larger is the number of coupled electronic states involved in the procesalj. But in non-linear spectra e-ph 
coupling can have much more peculiar effects with no counterpart in linear absorption. In conjugated materials 
e-ph coupling is responsible for th^firppearance in infrared and/or Raman spectra of vibrational bands with large 
(sometimes really huge) intensities! 1 -Hi 1 teII. Then there are vibrational states, to mean states within the ground state 
vibrational manifold, that have large transition dipole moments and/or transition (Raman) polarizabilities. These 
states, if properly included in sum over states (SOS) expressions of NLO susceptibilities, provide an important purely 
vibrational channel to NLO responses, that adds to the standard electronic channel. Vibrational channel contributions 
to resonant NLO processes have not been recognized so far: here we demonstrate they can strongly affect the spectral 
properties of 7r-conjugated systems. 

In the next Section we will concentrate on TPA spectra of push-pull chromophores described in terms of the simplest 
relevant model, the DA dimer with Holstein couplings. Within this model we are able to calculate numerically exact 
TPA spectra. We can then discuss the reliability of Condon and HT approximation schemes. We find sizable HT 
corrections to TPA spectra even for systems whose OPA spectra are well reproduced within the Condon approxima- 
tion. Even more importantly, we demonstrate large contributions of the vibrational channel to TPA spectra. These 
contributions become dominant in the cyanine limit. In the discussion Section we extend the treatment to other 
non-linear techniques and shortly address the behavior of different 7r-conjugated systems. 

II. TPA SPECTRA OF THE DA DIMER WITH HOLSTEIN COUPLING 

The simplest relevant model to describe (low-energy) spectral properties of push-pull chromophores accounts for 
just two-electronic states, \DA) and \D + A~), linearly coupled to a vibration, according to the following Hamiltonian 

(h = i)m 

H = z (l - a z ) - V2ta x - V2euj v Qp + ^{uj 2 ,Q 2 + P 2 ) (1) 

where Q and P are the coordinate and momentum for a harmonic vibration with frequency lo v ; o~ z and o~ x are the 
Pauli spin operators; 2zq is the energy gap at Q — 0; p = (1 — cr z )/2 is the ionicity operator measuring the amount 
of charge separation. The relaxation energy of \D + A~), e, measures the strength of the e-ph coupling. The above 
Hamiltonian assigns the two non- interacting (V2t = 0) states two harmonic potential energy surfaces (PES), as shown 
in Fig. la. Following Mulliken, we assume orthogonal basis states and neglect all dipole moment matrix elements but 
ji = (D + A~ \ jl\D + A~) , so that the dipole moment operator can be written as: /t = pojxQ- 

When the interaction is switched on (without loss of generality we set \J~2t = 1) the two states mix, and, in the 
adiabatic approximation, the ground and excited state PES, shown in Fig. lb, are easily obtained by diagonalizing 
the relevant two-dimensional matrix as a function of Q. As a consequence of e-ph coupling, i.e. of the relative 
displacement of the PES in Fig. la, the energy gap between the two basis states varies with Q, so that the amount 
of mixing between the two states, measured by p = (p), also depends on Q. Then, both the ground (\G)) and excited 
state (\E)) wavefunctions depend on Q, and analytical expressions for the Q-dependent properties (energies, dipole 
moments, etc.) are easily obtained. In particular the energies of both states show a non-parabolic dependence on 
Q, i.e. the two corresponding PES .are anharmonic (Fig. lb). This anharmonicity induced by e-ph coupling has 
large effects on static NLO responsesE3, where high order Q-derivatives of the ground state PES are involved. It also 
strongly afEects incoherent spectral measurements where the system tests large regions of the .ground and/or excited 
state PEScj. On the opposite, anharmonicity is not expected to show up in vertical processesES like OPA or TPA. 

In any case, the analytical expressions for the ground and excited state PES allow for an exact solution of the 
corresponding vibrational problemSj-A detailed description of the adiabatic solution of the Hamiltonian in Eq. (1) is 
deferred to a subsequent publicationEJ. Basically, chosen a harmonic phonon basis to describe the vibrational problem 
on either surfaces, the corresponding Hamiltonian matrices are easily written down by expanding the relevant electronic 
energies in powers of Q. Numerically exact vibrational eigenstates are then obtained provided that both the phonon 
basis and the Q-expansion are truncated at a sufficiently high order as to get convergence on the evaluated properties. 
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Similarly, the matrix elements of fi on the complete adiabatic basis, are readily obtained from the Q-expansion of 
the relevant electronic dipole moments (either pg(Q) = (G\fi\G), or pe{Q) = (E\p,\E) or pcr{Q) = (G\p,\E)) closed 
on the relevant exact anharmonic vibrational eigenstates. For the calculations presented in this paper, we chose 
to describe the vibrational problem in cither the ground or excited state PES on the basis of the eigenstates of a 
harmonic oscillator with frequency ui v , centered at the minimum of the ground state PES. With this choice we found 
that ~ 20 basis states and a similar number of terms in the Q-expansion of the electronic energies and/or dipole 
moments are more than enough to get convergence. We underline that the adiabatic approximation we adopt in the 
present treatment is very useful since it allows for a clear separation between states in the G- and E- manifolds, and 
then to distinguish between vibrational and vibronic excitations. The two kinds of states are completely mixed up in 
a truly non-adiabatic system. In any case, as far as the results presented in this paper are concerned, the reliability of 
the adiabatic approximatipH has been confirmed by the comparison with exact non-adiabatic results, obtained along 
the lines presented in Ref.ES. 

OPA spectra are calculated as wlmlx' 1 '^)], with the complex linear susceptibility, X ( w )> defined as: 

^H = E(^l/>(/IAI. 9 >(^ + ^) (2) 

where \g) — |G)|0) is the global ground state, |/) is an excited state, flf g = Uf g — iTf is the complex frequency of the 
/ .9 transition, defined in terms of the real transition frequency (u>f g ) and of the damping factor of the / state (Tf)cB. 
The summation in the above equation runs on all excitations, including vibrational excitations (in the G-manifold) as 
well as vibronic excitations (in the E- manifold) . The sum on vibrational excitations describes of course the infrared 
(IR) absorption spectrum, and basically does not contribute to the electronic spectrum, dominated by terms that can 
be written as (0\pct(Q)\v}(v\pct(Q)\0) /(^>ct,v — uj — iT), where |0) is the lowest exact anharmonic vibrational state, 
\v) is the v-th exact anharmonic vibrational state in the E'-manifold, and ll>ct,v is the |G)|0) — > transition 
frequency. 

In the following, we fix e = 1, and ui v — 0.1 as dime nsio n lfiss parameters (in \[2t units). In the specific case of 



push-pull chromophores, a typical value for \[2t is ~1 eVli3'E3c3, so that all the energy values can be approximately 
read in eV. We present results obtained for three different z values, to model chromophores with a very low polarity 
[zq = 0.95, p = 0.2), with intermediate polarity (zo = 0.66, p = 0.35) and just in the cyanine limit (zq = 0.50, 
p = 0.5). We notice that the properties of the system are symmetric with respect to the substitution p — > 1 — p. The 
upper panels in Fig. 2 report one-photon absorption (OPA) spectra: exact spectra (continuous lines) are compared 
to those obtained in the Condon (dotted lines) and HT (dashed lines) approximation schemes. Both Condon and HT 
schemes work .within the harmonic approximation, and we calculate Condon and HT spectra in the best harmonic 
approximation^ for the ground and excited state PES, i.e. by truncating the expansion of the corresponding potential 
energies, around the ground state equilibrium position, up to the second order in Q. In the Condon approximation, 
the expansion of pct(Q) is truncated at the zeroth order (constant) term, whereas in the HT approximation the 
linear term in the expansion of the numerator of x is also accounted for. As it turns out from Fig. 2, both Condon 
and HT approximations are fairly good for OPA spectra. The minor deviations of HT from Condon spectra are 
irrelevant in the fairly broad spectra usually observed for polar chromophores in solution. Quite predictably, the 
harmonic approximation works well to describe vertical electronic processesE3, and HT (and higher order) corrections 
to Condon spectra turn out small for allowed transitions. 

Single beam TPA spectra can be calculated as uj1oi[xtpa(~ uj, uj, — to)], where xtpa(— w; uj, uj, — u>) is the TPA 
component of the third order complex susceptibilityca: 

Xtpa (-uj;uj,uj,-uj) (x ^2 (fflAl TO )( m lMl/)(/lMl«X«lAl.9> x 

i i 

+ 



(Q mg - uj)(Qf g - 2uj){VL ng - uj) (fi; ng + uj)(Q fg - 2uj)(Q ng - uj) 
1 1 



(n* mg + uj) (n* fg + 2u) (n ng - uj) (n* mg + uj)(n* fg + 2uj) (n* ng + uj) 



(3) 



where the triple sum runs on all excited states. Of course, if we are interested in TPA processes leading to the 
E- manifold, (i.e. processes measured for uj ~ ujct/2), then / only counts vibronic excitations. But m and n run on 
both vibronic and vibrational excitations. ____ 

In all available studies of vibrational contributions to TPA spectraD&Efli3, all summations in Eq. (||) have been limited 
to vibronic excitations only, i.e. only the electronic channel to TPA has been accounted for. In this approximation, 
Xt pa reduces to: 
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E (0\»ct(Q)W)(u\he 

u,v,w 

1 + 1 

(&CT,u - u)(QcT,v - 2u)(QcT,w - <Jj) (^CT.u + W )(^C7> ~ 2w)(fic2> - to) 



("cr,™ + w) (fi^ + 2w) (0 C r,. - w) (ft CT> „ + co) + 2w) (O*,^ + w) _ 

where tt,i> and w run on vibrational states in the iJ-manifold, and p — (0|^g(Q)|0) is the (Q-independent) ground 
state dipole moment. In the middle panels of Fig. 2 we report the electronic channel contribution to TPA spectra, 
calculated for the exact adiabatic eigenstates of the Hamiltonian in Eq. (1) (continuous lines). In the same figure, 
dotted and dashed lines refer to spectra obtained within the Condon and HT approximation schemes. 

Even accounting for just the electronic channel to TPA absorption, we can already stress the importance of a proper 
modeling of e-ph coupling in TPA spectra. If e-ph coupling is not accounted for, both OPA and TPA spectra have a 
single Lorentzian shape, and are exactly superimposed if the TPA frequency scale is expanded by a factor of 2 with 
respect to the OPA scale. At the Condon level, the situation is not very different: e-ph coupling originates structured 
bands in both OPA and TPA spectra, but with similar shape, so that the two spectra are essentially superimposable 
as in the absence of coupling. An obvious deviation is observed at p = 0.5, where, in the harmonic approximation, 
the TPA spectrum exactly vanishes due to the factor pe — f-i, proportional to 1 — 2p in the dimer model. At variance 
with OPA spectra, however, HT corrections are very important in TPA, at least for not too neutral molecules. For 
p = 0.35, the HT spectrum in the central panel of Fig. 2 has a qualitatively different shape from the Condon spectrum, 
with the 0-1 line acquiring intensity at the expense of the 0-0 line. Apart from the pathological p = 0.5 case, the HT 
approach offers a good approximation to the electronic channel contribution to TPA, confirming the reliability of the 
harmonic approximation to describe vertical electronic processesE3. 

Up to now we have only considered the electronic channel contribution to TPA process, where e-ph coupling mainly 
contributes in defining vibronically structured band-profiles. We already got a fairly impressive result: the Condon 
approximation to TPA spectra is in general inadequate, even in cases where it works well for OPA spectra. However 
in TPA spectra, and more generally in non-linear spectra, vibrations have a much more important role than simply 
structuring the observed band-shapes. They in fact open a vibrational channel to non-linear processes. The lower 
panels in Fig. 2 show the total TPA spectrum, obtained from the imaginary part of xtpa in Eq. (^), by allowing m 
and n indices to run on both vibrational and vibronic excitations. The comparison between continuous lines in the 
middle and lowest panels in Fig. 2 clearly demonstrates the importance of the vibrational channel to TPA process, 
with sizable effects on the absorption frequency, on the band-shape as well as on the total intensity. Overall, the 
vibrational channel increases the intensity of higher vibronic replicas at the expense of the lowest replicas, with effects 
that are most important at p ~ 0.4. This implies a blue-shift of the TPA spectrum if compared with the rescaled 
OPA spectrum. In the extreme case of p = 0.5, the blue-shift amounts to a vibrational spacing. Intensity effects are 
negligible for p ~ 0.2, but at p — 0.35 the vibrational channel amplifies TPA response by a 1.2 factor, that becomes 
a factor of 4 in the p = 0.5 case. 



III. DISCUSSION 



Vibrational channel contributions to TPA process arise from terms in Eq. (0) where m or n, or both m and n run on 
vibrational excitations, so that in each term contributing to the vibrational channel at least a factor like (0\pg{Q)\v ) 
appears, where v is one of the vibrational states in the G-manifold. The square of this term measures the absosption 
intensity of the v-th vibrational state. Since the fundamental vibration has by far the largest IR intensityE3, the 
vibrational channel to TPA is dominated by it. A first approximation to Xtpa ^ s then: 



vib 

Xtpa 



MR(MtiCT(Q)\v)(v\f*E(Q) ~ p\w)(w\pct(Q)\0) 
{U) v - Lj){VtcT,v - 2lu)(Qct,w - w) 



E 



MJgCjj MCT {Q)\v) {v\VCT (Q) 1 1) HIB, 

(lo v - lu)(Q C t,v - 2lu)(uj v - to) 



(•5) 



where pm = (0\pq(Q)\1) is the dipole moment of the fundamental vibrational transition. The relative importance of 



Xtpa vs Xtpa g° es as - A*)) Oct 



;|) or as its square, depending if the first or the second term 
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in the above equation is concerned. Since (Xir oc [p(l — p)] 3 / 2 @ maximizes at p = 0.5, whereas (1 — 2p) 

vanishes there, the relative importance of Xtpa increases fast towards the cyanine limit. On the other hand, the 
frequency of vibrational modes, u> v ~ 0.1 — 0.2 eV, is not far from the laser frequency to ~ u>ct/2 ~ 0.5 — 1.5 eV, so 
that the ratio {{ujct — &)/(k>v ~~ <^)) enhances the vibrational contribution over the electronic one by a factor ranging 
up to 3. 

The solution of the two-state model in the adiabatic approximation allows for reliable estimates of the relative 
magnitude of electronic and vibrational contributions to TPA spectra and even more importantly offers a natural 
understanding of the physical origin of the different contributions. An obvious drawback of the two-state model is 
however the impossibility to account for the role of highly excited states. Of course, extending the sums in Eq. (||) to 
excited states (E') higher than the lowest one (E) opens additional channels to the TPA process. These .E'-channel 
contributions to TPA have m or n or both m and n in Eq. (3) running over vibronic states pertaining to the E'- 
manifolds. By closing over the vibrational states of these high-energy levels, one can roughly estimate the magnitude 
of S'-channel contributions in terms of the G — ► E' and E — > E' transition dipole moments (beside of G — > E) and of 
the corresponding transition energies. At the best of our knowledge, experimental estimates of the E — > E' transition 
dipole moments are only available for coumarin 153 dispersed in CH3CNED. For this molecule, excited state transient 
absorption measurements yield a E — > E' transition dipole moment of ~ 10 D. Transition dipole moments of similar 
magnitude are also evaluated, from OPA-spectra, for the two transitions G — * E and G — > E' , occurring at energies 
of ~ 3.0jsy and ~ 5.6 eV, respectivelyeZlEa. The joined analysis of electronic and vibrational spectra of the same 
moleculeEEl yields p*c* 0.23, e ~ 1 eV and ui v ~ 0.2 eV, in fairly good agreement with previous estimates on a similar 
dye (coumarin 102)113. Based on these data, the vibrational and i?'-channel contributions to TPA spectra turn out of 
the same order of magnitude. In view of the fairly low polarity of coumarin 153 and of its fairly high loqt frequency, 
we believe that for most push-pull chromophores the vibrational channel contribution to TPA spectra is larger than 
£"-channel contributions. In fact, preliminary ZINDO calculations on a few chromophores confirm that, at least for 
chromophores with intermediate polarity (p ~ 0.3 — 0.7), the vibrational contributions to TPA spectra are important 
and cannot be disregarded with respect to .E'-channel contributions. Refined models for NLO responses accounting 
for plenty of excited states are bound to fail if no proper account is given for the contributions of vibrational states. 

Besides TPA spectra, other NLO responses can be affected by the vibrational channel. Fig. 3 shows the amplitude 
of the third harmonic generation (THG) signal in the region uj ~ wct/3, for three different p values. Continuous lines 
refer to the total response, obtained by allowing the triple sums in xthg = X ( — 3cj; u, oj, lo) to run on both vibronic 
and vibrational excitations. Dashed lines instead account for the electronic channel to THG, obtained by summing 
only on vibronic states. The vibrational channel strongly affects THG response, with largest effects at p ~ 0.3, where 
the response is amplified by a factor ~ 2 in the shown spectral region. Calculations are performed using the complete 
expression for xthg-, but, in the portion of the spectrum shown in the figure, the terms with a resonating denominator 
are clearly dominant, so that we can approximate xthg as follows: 



Xthg(Su; u, w,u) ~ 

9\fr\f)(f\l*\m)(m\jl\n)(n\p,\g) \ - (s|A|/)(/|A|s)(3l£l™)( m IA|s) 



E \y\p\j 1 \j v'"-imi'v K ii myi _ 



, ni n ( n fg ~ 3w ) i^mg - 2w) (O ng -w) jr£ {Q fg - 3w) (Q fg - u) (Cl mg - w) 



(6) 



where / runs on vibronic excitations, just to guarantee an almost resonant denominator, luct ~ 3w, and m,n run 
on both vibrational and vibronic excitations. The vibrational channel contribution to THG is more complex than 
for TPA. A few terms of Xthg-> obtained by allowing n or m or both m and n in the above equation to run on 
vibrational excitations, contain at the numerator fim or its square. A term however has a different numerator: 
J2 V w(^\f 1 CT(Q)\v}(v\pcT(Q)\i}(Mt J 'CT{Q)\w)(w\ficT{Q)\0), that is related to the non-resonant Raman intensity of 
the fundamental vibration. It is difficult to disentangle the contributions of the various terms to the observed signal, 
but the importance of Raman modes is suggested by the observation of the maximum vibrational contribution to 
THG just at p ~ 0.3, where the non-resonant Raman intensity maximizes for the DA dimer modelcJ. 

Having rationalized the vibrational channel contributions to TPA and THG spectra of the DA dimer model in 
terms of IR and non-resonant Raman intensity of coupled vibrations, we can shortly extend our discussion to different 
systems. In conjugated polymers the coupled vibrations are not IR active, whereas they have huge (non-resonant) 
Raman intensitiesErtj. Then we expect no contribution from the vibrational channel in TPA spectra of pristine 
polymers, whereas the vibrational contributions to THG spectra can be important. A detailed description of NLO 
properties of polymers is difficult, in view of their [fairly crowded excitation spectra. Preliminary results on a simple 
relevant model, the self-dimer with Peierls couplingEj, confirm large vibrational channel contributions to THG spectra, 
and negligible contributions to TPA spectra. Of course large vibrational contributions to TPA spectra can be expected 
for doped or photocjajitpd polymer samples, based on the very large IR intensities associated with charged solitonic 



or polaronic defectaij'cil. 
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IV. CONCLUSIONS 



In conjugated molecules and polymers e-ph coupling is very effective_.lt shows up with largely anharmonic PES, 
that, in conjugated polymers, often show more than a single minimumtj to signal the Peierls instability. A second 
fairly ubiquitous consequence of e-ph coupling is the large nuclear, displacement upon electronic excitation, with 
well apparent effects in absorption and resonant Raman spectraB— ]■— lr 6 H 29 l. In systems with large e-ph coupling, 
static electric fields can draw nuclei out of equilibrium, and this is the basic reason why e-ph coupling is so effective 
in amplifying static NLO responses via anharmonic contributions— I. Fields at optical frequencies are too fast to 
displace nuclei, but in incoherent measurements nuclei are allowed to relax for a while following the ap p l ic ation of the 
preparation pulses: once more e-ph coupling is important in defining the shape of the involved PES— Electron- 
phonon coupling largely affects vibrational spectra, with a softening of the relevant modes and a large increase of 
infrared and/or Raman intensities. Therefore in SOS expressions for the susceptibilities, vibrational states can give 
large contributions. For static susceptibilities this contribution indeed accounts for the anharmonicity of the ground 
state PES_H_. 

As far as resonant processes are concerned, the importance of the vibrational contributions has not been fully 
appreciated so far. Here we have analyzed TPA spectra of the DA dimcr model — the simplest model relevant to 
push-pull chromophores — and proved that the vibrational contribution is very important in defining the shape, 
frequency and intensity of the signal. Including vibrational degrees of freedom has two major effects in TPA spectra, 
and more generally in resonant NLO processes. TheJjkstpJjairly obvious, effect is the appearance of a vibronic structure 
in TPA spectra. This effect was already recognize<fii§Bt_; here we prove the validity of the harmonic approximation 
in this context, but the inadequacy of the Condon approximation. Specifically, HT corrections to TPA spectra are 
important even in cases when OPA spectra are well reproduced within the Condon approximation. The second effect 
of vibrations is more interesting and has not been recognized so far. Vibrational states in the ground state manifold 
in fact contribute a new and important channel to TPA process. The vibrational channel contribution shows up 
with a blue-shift of the TPA spectrum with respect to the rescaled OPA spectrum and with a large distortion of the 
TPA band-shape. These contributions are related to simple spectroscopic observables as inferred from vibrational 
and linear absorption spectra. This immediately leads to safe predictions about other systems. As a matter of fact, 
vibrational channel contributions to TPA spectra are related to the IR intensity of coupled modes, so that they are 
expected to be negligible in any polymer (but possibly dominant in doped and/or photoexcited polymers). THG 
spectra instead have contributions from Raman modes, so that large effects of e-ph coupling are expected in THG 
spectra of conjugated polymers. 

The behavior of conjugated materials can only be understood by properly accounting for e-ph coupling, that affects 
in different ways and for different reasons the spectral properties of these systems. Common beliefs based on the 
behavior of other systems where e-ph coupling is not so effective, and/or extrapolated from our knowledge of their 
linear spectroscopic behavior need to be carefully reconsidered. A joint experimental and theoretical effort is needed 
to extract the correct physics out of these fascinating materials. 
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FIG. 1. Potential energy surfaces sketched along the vibrational coordinate Q for (a) the basis states {\/2t = 0), and (b) 
the exact eigenstates (V^t — 1). Both panels refer to zq = 0.95 and e = 1. 
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FIG. 2. OPA and TPA spectra calculated for e = 1 and oj v — 0.1 in units of \/2t. Left, middle and right panels refer to 
p = 0.2, 0.35, 0.5, respectively. The vibronic and vibrational damping factors are fixed to 0.05 and 0.005, respectively. Upper 
panels report exact OPA spectra as continuous lines. Dashed and dotted lines correspond to spectra obtained in the HT and 
Condon approximation, respectively. Middle panels report the electronic channel contribution to TPA spectra; continuous, 
dashed and dotted lines refer to exact spectra and to spectra obtained in the HT and Condon approximation, respectively. 
Bottom panels show the total TPA spectrum including the vibrational channel contribution. The a:-axis scale at the top and 
bottom of the figure report the energy of the OPA and TPA spectra, respectively. j/-axes (arbitrary units) preserve the relative 
intensities. 
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FIG. 3. THG spectra calculated for the same parameters as in Fig. 1. Upper, middle and bottom panels refer to p = 0.2, 
0.3, 0.5, respectively. THG spectra including the vibrational channel contribution (continuous lines) are compared with the 
purely electronic THG spectra (dashed lines). 
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